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Abstract 

An analytic calculation of the photon mass gap M of compact U(l)2+i 
in the Hamiltonian formalism is performed utilizing the first four Hamil- 
tonian moments with respect to a one-plaquette mean field state in the 
plaquette expansion method. Scaling of M is clearly evident at and be- 
yond the transition from strong to weak coupling. The scaling behaviour 
agrees well with the range of results from numerical calculations. 

1 Introduction 

In QCD, the most fundamental quantities one wishes to calculate are the states 
of the mass spectrum of the theory. In lattice gauge theory, such non-perturbative 
calculations are carried out numerically in the path integral representation to 
some level of approximation by evaluating the correlation functions by Monte- 
Carlo simulation. Typically, the computational resources employed in this ap- 
proach are enormous as one must work with large lattices and spacings which 
are fine enough to detect the continuum behaviour of the theory in question. 
Hence, non-perturbative analytic results are of considerable interest. In this 
context, the Hamiltonian formalism is appealing. One often works directly on 
infinite lattices, fermions are simpler to incorporate (no fermion determinant) 
and 3D as opposed to 4D lattices (time is continuous). The lack of systematic 
and applicable methodology has prevented these features from being exploited 
or investigated - hence, the importance of new methodologies. Here we use a 
novel approach - the plaquette expansion - based on a large volume expansion 
of Lanczos tridiagonalization. 
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2 Compact U(l) in 2+1 Dimensions 



Compact U(l)2+i has become the testing ground for various lattice Hamiltonian 
procedures as it exhibits scahng behaviour similar to the more complex and 
physically interesting non-abelian lattice gauge theories in 3 + 1 dimensions. 
The Kogut-Susskind Hamiltonian is |^ : 

^-yE^' + ^E[2-(^^P + ^;)] - (1) 
I ^ p 

where g is the dimensionless coupling constant. The strong-coupling limit is 
defined by g — > oo and the weak-coupling limit hy g 0. The electric field 
operator, Ei, and link operator, Ui obey the commutator relations [Ei, Ui] = Ui 
and [El, Uj] = — fj} and the plaquette operator (jp acts on the links around the 
smallest closed (Wilson) loop or square on the lattice JJp — UiU2Ulu\- 

One of the quantities of interest in this model is the anti-symmetric or photon 
mass gap, M ( = M^+- ), which is given by the difference in energies between 
the lowest state in the 0"' sector and the vacuum. The scaling behaviour for 
the mass gap is expected to be M'^a? = /?exp{— fco/9 -I- ki\, where and 
fci are constants, a is the lattice spacing and (3 — l/g^. The scaling parameters 
fco and fci are not known exactly but have been computed by various numerical 
techniques: fco ^ 4.1 — 6.345 and ki ^ 4.369 — 6.27 (See p| for summary). 



3 Plaquette Expansion 

Beginning with a trial state \tpi) which has the desired symmetries of the state 
of interest, the Lanczos recurrence generates a basis 



f3n- 



[{H - a„_i) \tpn-l) - l3n-2\lp7i-2)] 



where a„ = {'ipn\H\ipn) and /3„ — {iJjn+i\H\'ipn) are the matrix elements of the 
Hamiltonian in tri-diagonal form. 

The matrix elements are able to be written in terms of Hamiltonian mo- 
ments (i?") = {ipi\H"\il^i) . The connected part of the Hamiltonian moment is 
proportional to the volume of the system {H")c — SnN. Therefore, one may 
re-express the matrix elements in terms of the connected vacuum coefficients 
e„. Extensivity of the problem leads to the following plaquette expansions in 
1/iV i: 



On 

N 



ei 



Pi 

iV2 



= se2 + -s 



£3 
6264 



+ 0{S^) EE a(s), 



(2) 



(3) 



where s = n/N. In the bulk limit (n, N — > oo), keeping s fixed, one may perform 
the exact diagonalization of the Lanczos tri-diagonal matrix for the ground state 
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energy density analytically Q: eg = inf (a(s) — 2/3(s)) . For example, for 4th 
order moments we have: 



g2 

eo[4] = ei + 2-^ 

6264 — 63 



3ei — 26264 — 63 



(4) 



{ s) 

For excited states the Hamiltonian moments have the form {H^)c = CnN + 
Thus, in a similar fashion to £0, one can derive expressions for approx- 
imants, Afg [r], to the mass gap. Again using up to 4th order moments we 
have |6j: 
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Mf )[4] - ™« ^ ^2"' K(e2, 63, 64), (5) 
where the vacuum moment functions, Fn, are simple algebraic forms given in Q. 

4 One Plaquette Mean-Field State 

Initially, lattice Hamiltonian calculations employed the simplest gauge invari- 
ant trial state — the strong-coupling vacuum |0) (defined as the state satisfying 
Ei\Q) = for every link). This is the perturbative starting point for series 
calculations, which are then extrapolated to weak coupling. Although the cal- 
culation of moments with respect to this state can be carried out to relatively 
high order, typically, one finds that this state is simply inadequate to explore 
the non-perturbative weak-coupling regime of the theory. In the case of the 
plaquette expansion, a window of scaling in terms of the expansion order was 
evident, but higher order results became problematic 0. 

Here we utilized a one plaquette exponential trial state, that has the correct 
strong and weak coupling behaviour for the vacuum energy density - the so- 
called mean-field state - defined in the so-called U representation as: 

V^iif/J^expMviC/p + C/Jjl (6) 



where A is determined variationally from {H) . As the state is constructed from 
plaquette variables. Up, it is automatically gauge invariant. One expects the 
true ground state of the system to have contributions from all possible Wilson 
loops - an important feature of the plaquette expansion method is that, by 
its very nature, larger sized Wilson loops (larger clusters) are systematically 
introduced. 

Using a diagrammatic method (details to be presented elsewhere) we were 
able to derive the first four one-plaquette connected moments for both the 
ground state and anti-symmetric state. The expressions become rather long, 
so for brevity these are illustrated in Figures |l| and |2[ along the curve A(/3). 
These moments immediately give analytic approximates for vacuum energy and 
the mass gap. 



3 



6.0 



5.0 



4.0 



3.0 



2.0 



1.0 



0- 



1 1 1 


1 

•■•*■. 

■ * * 


1 1 1 


• 


• 


, — 

'■•..^3 _ 


• 










eg - 






1.0 3.0 3.0 4.0 

Figure 1: Ground state moments, ei 64. 
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Figure 2: Anti-symmetric moments, mi — > 7724. 



5 Results and Discussion 



The plaquette expansion results are presented in Figure ^. Clearly scaling is 
evident for inverse-coupling values /3 = 0.7 to 1.5, passing the transition point 
at /3 = 1.0. The scaling form we find is given by = /3exp (-4.34/3 -I- 5.34), 
and agrees well with other estimates (summarized in [D). 
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Figure 3: Mass gap results applying one, three and four one-plaquette Hamilto- 
nian moments. Scaling is evident for the four order moments curve from /? — 0.7 
to 1.5. 

Although we have presented results for only the first four moments (for which 
the calculation is entirely analytic) it should in principle be possible to derive 
higher moments with respect to this state. Indeed we have preliminary results 
for the ground state energy density to sixth order giving high accuracy over a 
large range of couplings. 

For the future application of this method to QCD3_,.]^, the main obstacle 
to overcome is the calculation of the moments in 3-1-1 dimensions. The appeal 
of 2+1 systems as toy models is that the transformation from link variables 
to plaquette variables is trivial and makes the integrations tractable - hence 
the analytic work reported here. In 3-1-1 dimensions it is well known that the 
transformation cannot be carried out in closed form due to the appearance of 
Bianchi identities. However, the calculation of moments in 3-1-1 dimensions by 
Monte-Carlo is actually a relatively small scale numerical exercise. Because we 
are dealing with a cluster expansion (we need only connected moments) the lat- 
tices required are fairly modest in extent and, by definition, only 3 dimensional. 
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The integrands become complicated as larger correlations arc included, which 
means that the statistics must be very good, however, preliminary calculations 
have demonstrated that the required precision is possible to achieve without 
supercomputing resources. 
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